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On the Redshift of the Very High Energy Blazar 3C 66A 

A. Furniss^, M. Fumagalli^'^''', C. Danforth^, D. A. Williams^, J. X. Prochaska^ 

ABSTRACT 

As a bright gamma-ray source, 3C 66A is of great interest to the high-energy 
astrophysics community, having a potential for placing cosmological constraints 
on models for the extragalactic background light (EBL) and the processes which 
contribute to this photon field. No firm spectroscopic redshift measurement has 
been possible for this blazar due to a lack of intrinsic emission and absorption fea- 
tures in optical spectra. We present new far-ultraviolet spectra from the Hubble 
Space Telescope/Cosmic Origins Spectrograph (HST/COS) of the BL Lac object 
3C66A covering the wavelength range 1132 to 1800 A . The data show a smooth 
continuum with intergalactic medium absorption features which can be used to 
place a firm lower limit on the blazar redshift of z > 0.3347. An upper limit 
is set by statistically treating the non-detection of additional absorbers beyond 
z = 0.3347, indicating a redshift of less than 0.41 at 99% confidence and ruling 
out z >0.444 at 99.9% confidence. We conclude by showing how the redshift 
limits derived from the COS spectra remove the potential for this gamma-ray 
emitting blazar to place an upper limit on the flux of the EBL using high energy 
data from a flare in October of 2009. 

Subject headings: BL Lacertae objects: individual: 3C 66A, galaxies: active, 
intergalactic medium, ultraviolet: general, cosmology: diffuse radiation 
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Direct measurements of the extragalactic background light (EBL) are difficult due to 
stron g foreground sources in our solar system (Zodiacal light) and the Galaxy flHauser fc Dwek 
200X1 ). If a direct measurement were possible, it would only reflect the current integrated 
state, leaving still the model-dependent task of extracting the time evolution. These dif- 
ficulties have been overcome through the use of extragalactic ver y-high-energy (VHE, E> 
100 GeV) garn r na ray s from blazars with k nown redshifts f e. g. Aharonian et al.l (2006 ): 
Mazin fc Rau3 J2OO7I ): klbert et al.l fcoosh : Icilmore et all J2OO9I ): lOrr et al.l fcoilh ). the 
most commonly detected type of VHE extragalactic source. A blazar is a type of active 
galactic nucleus (AGN) that has a jet pointed towards the observer, and exhibits a highly 
polarized broadband spectrum from beamed, non-thermal emission processes. 

The energy-dependent absorption of gamma rays by the EBL softens the intrinsic VHE 
gamma-ray spectra emitted by extragalactic objects. The details of the absorption de- 
pend on the shape of the EBL spectral energy distribution (SED) in the near-IR to optical 
band. Additionally, the total power and the shape of the SED of the EBL is shown t o var y 
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the gamma-ray absorption, an accurate redshift of the VHE extragalactic target is required. 



Approximately one-third of the current VHE extragalactic catalogj is made up by 
blazars at unknown or poorly-constrained redshift. 3C66A is one of these blazars, with 
an un certain spectroscopic redshift based on possible corrobor ating measurements of single 
lines ( iMiller et al.lll978l : iKinney et al.lll991t iBramel et al.ll2005[ ). Despite multiple attempts, 
in particular two high signal-to- noise measurements using Keck I (shown in Figure 1), no 
solid spectroscopic measurement based on the detection of multiple lines from the host galaxy 
has been possible. The lack of spectral features is not surprising given that 3C 66A is a BL 
Lac-type active galactic nucleus that, by definition, displays weak or no lines. 

To overcome the inherent featureless characteristic of the 3C 66A optical spectrum and 
enable deabsorption of the VHE spectrum with reliable redshift information, we have deter- 
mined limits on the redshift of the blazar through the observation and statistical analysis of 
far UV (FUV) absorption by the low z intergalactic medium ( IGM). This method, alread y 
applied to the VHE blazars PG 1553-M13 and S5 0716+714 dDanforth et^ boiol . l2013h . 
sets a redshift lower limit using absorption lines caused by the intervening IGM. Further, 
given the expected distribution of IGM absorbers as a function of redshift, one can model 
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40 any lack of absorption lines at longer wavelengths to statistically infer an upper limit on the 

41 blazar redshift. 



42 2. Observations and Spectral Analysis 

43 3C 66A was observed with the Cosmic Origins Spectrograph (COS) during two epochs as 

44 part of two different programs. The blazar was observed for five HST orbits on 1 November 

45 2012, with the medium resolution G130M grating COS/G130M (1135 < A < 1450 A, 15.3 

46 ksec) as part of program 12621 (PI: Stocke). Three more orbits were devoted to observations 

47 with the COS/G160M (1400 < A < 1795 A, 7.2 ksec) grating under program 12863 (PI: 

48 Furniss) on 8 November 2012. The calibrated, one-dimensional spectra for each exposure 

49 were obtained from the Mikulski Archive for Space Telescopes (MAST). 

50 The G130M data show a flux mis-match between the short and long-wavelength segment 

51 of each exposure and a ~ 8% correction is applied to each before coaddition to bring them 

52 into the expected smooth continuum. The G160M data are considerably noisier; no flux 

53 discrepancy was observed and no correction was undertaken. The co rrected exposures wer e 



then coadded with the standard IDL procedures described in detail by lDanforth et al.l f l2010l ). 
This procedure includes an automatic scaling of the exposures taken during different epochs. 
The continuum flux level appears to have varied by < 10% during the week between observing 
epochs, well within the current flux calibration uncertainty. 

The combined spectrum continuously covers the wavelength interval 1132 — 1800 A , and 
shows the expected smooth continuum and narrow absorption features. The data quality 
varies over the spectral range due to the different sensitivities and exposure times in the 
two gratings. The mean signal to noise per pixel in the continuum is ~ 10 (~ 5) with 
nominal dispersions of 9.97mA /pixel (12.23A /pixel) in the G130M (G160M) portion of the 
spectrum. Signal to noise value s per s even-pixel resolution element are approximately twice 



64 these values (see Keenev et al.l f|20 



Ghavamian et al.l (120091 ) and 



Kriss 



1) ) . F or additional details on the COS instrument, see 



torn . 



Detailed analysis of these data and of the intervening absorption line systems will be 
presented in Danforth et al. (in prep.). In this paper, we exclusively focus on the spectral 
features that are useful to constrain the unknown redshift of the blazar (2;biazar)- The goal 
of the following analysis is to use absorption lines that arise from gas clouds in the IGM 
to establish a firm lower limit on the distance to 3C 66A and to set an upper limit for the 
blazar redshift based on a statistical argument. 

A visual inspection of the spectrum reveals the presence of multiple absorption systems 
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for which both Lyman— a and Lyman—/? (Lya and Ly/3) hnes are detected. Among those, 
we identify three clouds at Zabs ~ 0.3283, 0.3333, and 0.3347 (see Figure [2]). All other lines 
detected at >4 sigma significance redward of these three Lya systems are identified as Milky 
Way absorption (see Figure [3]). Thus, because of the presence of a system a.t z = 0.3347, 
we set a secure redshift lower limit of 3C66A at ^biazar > z\\ = 0.3347. We also search the 
spectrum for O VI (AA1031, 1037) doublets that, owing to their bluer rest-frame wavelengths, 
could yield a more stringent redshift lower limit than the one set by absorption in the Lyman 
series. However, we do not find any instances of absorption beyond z ~ 0.33. 

Next, we can exploit the lack of absorption beyond zn = 0.3347 to set an upper limit 
Zui to the blazar redshift following a statistical argument. The frequency of absorption lines 
arising from the Lyman forest in the local uni yerse has been measured along sightlines t o 



extragalactic sources by different authors (e.g. iPenton et al.ll2004l : iDanforth fc Shullll2008l ). 
It is common to express this quantity with the function dN{W > Wo)/dz which describes 
the average number of absorption lines with rest-frame equivalent width in excess to Wq per 
unit redshift. We can therefore estimate the number of lines we expect to detect between 
z\i and z^i, given the rest-frame limiting equivalent width of the COS spectrum. By 
comparing the predicted number of absorption lines in a given redshift interval with the lack 
of detection beyond zw = 0.3347, we obtain redshift upper limit. 

First, we generate 1000 mock spectra in the observed wavelength range 1215 — 1800 
A by drawing Lyman forest lines from a distribution as a function of redshift such that 
the number of lines satisfies the observed dN{W > Wo)/dz. In this analysis, we assume 
no evolution in Lyman- q forest line incidence and adopt the frequency distribution from 
Danforth fc Shulll (120081) . although a similar result is obtained if we adopted the distribution 
from lPenton et al.l (120041 ). Next, we as sign to each line a Dopp ler parameter drawn from the 
observed distribution in the local IGM (IDanforth &: Shullll2008l ). During this step, we assume 
that the Doppler parameter is not correlated with the equivalent width of the line. Given 
a line equivalent width, its redshift, and a Doppler parameter, we compute the observed 
limiting equivalent width W^ra (at 5a) using the formalism developed for COS spectra by 
Keeney et al.l (120121 ) and we record only those lines which would be detected in the observed 
COS spectrum. Note that this procedure naturally accounts for "shadowing" due to Milky 
Way absorption lines. 

The top panel of Figure H] shows the number of intervening absorption lines detected in 
1000 mock spectra within the redshift interval 0.335 ^ z < 0.444. According to this Figure, 
we should expect to detect ~ 5 or more lines if 3C66A lies at ^ui > 0.444, and, although 
realizations with no lines are possible, they are extremely rare (< 1% of the total trials). 
Under the simplistic assumption that the number of absorption lines is not correlated in 
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velocity space, the mock realizations shown in the top panel of Figure 4 follow a Poisson 
distribution. Therefore, we can adopt Poisson statistics to express the probability of finding 
no detected lines between zw and z^i, given a typical number of Lyman forest lines in that 
redshift interval iV(zu < z < z^i). 

As shown in the bottom panel of Figure HJ the expected number of absorption lines 
increases proportionally to the redshift interval Az = z^i — zn and the probability of finding no 
absorption lines P{N = 0) exponentially decreases with redshift. At ^ui ~ 0.41, P{N = 0) ~ 
0.01, and therefore we conclude that 3C66A is likely to lie between 0.3347 < -Zbiazar ^ 0.41. 
We can further rule out ^biazar ^ 0.444 based on the fact that P{N = 0) ~ 0.001 for 
Zui ~ 0.444. We note that consistent probabilities can be recovered directly from the Monte 
Carlo simulations, without explicitly using Poisson statistics. However, it should be noted 
that our Monte Carlo simulations do not include correlated absorption systems in the Lya 
forest. Further, this calculation does not account for mechanisms that could enhance (e.g. 
galaxy clustering) or suppress (photoionization along the line of sight) the incidence of Lya 
lines in proximity to a blazar compared to the mean value observed in the IGM, although 
there is no evidence for highly ionized gas ( i.e. N V absorption) at 2 ~ 0.335. With 
proximity effects included, the predicted limits are subject to ~ 1000 km s~^ uncertainty 
(i.e. ~ 0.003 in redshift space). Nota bly, there have beeri previous suggestions tha t 3C66A 
is a member of a cluster at z ~ 0.37 f lButcher et al.lll976l : IWurtz et al.lll993l . 119971 ). 



The limits placing the re dshift betwe e n 0.3 347 a nd 0.41 disfavor the past tentative 
measurements oi z = 0.444 by iMiller et al.l ( 119781 ) and iKinney et al.l f ll99ll ). both of which 
were based on the measurement of single, weak lines. The limits derived from the COS 
observations are, h owever, in good agreement with other past estimates of the blazar distance. 
Finke et al.l (120081 ) set a lower limit of z > 0.096, an estimation based on the expected 
equivalent widths of absorption features in the blazar host galaxy, while a distance estimate 
of 2 ^ 0.321, noting a large error, was formed based on the assumption that host galaxies 
of a blazars could be taken as standard candles. An estimate for the blazar redshift of z = 



0.34 ± 0.05 was found by iPrandini et al.l (120101 ) . who extracted the approximate redshift by 
correcting the TeV spectrum of the blazar for EBL absorption t o match the index measured 
by the Fermi Large Area Telescope (LAT; lAtwood et al.l ( 20091 )). most sensitive to gamma 
rays between 300 MeV and ~100 GeV which are largely unaffected by the EBL. The redshift 
limits for 3C66A are also in good agreement with a recent EBL model-independent study 
of the gamma-ray horizon, as determined by syn chrotron self-Compton modeling of VHE 
blazar broadband spectra ( jPominguez et al.ll2013l ). 
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3. Absorption of Very High Energy Gamma-rays from 3C 66A 
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The energy- and redshift-dependent absorption of gamma rays by the EBL can be 
estimated using model-specific optical depths, t{E,z), where the intrinsic flux (Fint) can be 
estimated by the observed flux (Fohs) using the relation F^nt ~ Fobs x e^^^'^-*. The intrinsic 
index of a blazar can be used to estimate the spectral properties of the EBL under the 
physically motivated assumption that the intrinsic spectrum of a source undergoing Fermi 
shock acceleration, characterized by the power-law dN/dE cx , cannot be harder than 
r = 1.5. If the intrinsic VHE spectrum is signiflcantly harder than the r=1.5 limit, it can 
be argued that the gamma-ray opacity of the EBL model which was used for deabsorption is 
too high. The index limit of L5 is derived from the standard leptonic and hadronic emission 
scenarios used to describe blazar non-thermal emission. This limit is als o in agreement wit h 
the hardest gamma-ray index reported by the Fermi LAT for a blazar flNolan et al.ll2012l ). 
The indices for sources derived from photons with energies of less than 100 GeV are not 
signiflcantly affected by EBL absorption and so reflect the intrinsically emitted spectra of 
blazars in the high energy gamma-ray band. 

Under the assumption that blazars do not harden with increasing energy, EBL flux con- 
straints are also possible by comparing deabsorbed VHE spectra to the extrapolations based 
on the LAT-measured spectral indices. Using this method, the Fermi and VERITAS indices 
measured during a state o f elevated flux frorn 3C 66A in October of 2009 (F = l.SiO.l^jaj and 
4.1 ± O.Gstat, respectively; lAbdo et al.l (1201 if )) allow the investigation of possible constraints 
on the EBL density, pending a reliable distance measurement. Previously, the deabsorption 
of the VHE spectrum o f 3C 66 A has been c ompleted with the unce rt ain spectroscopic re d- 
shift of z = 0.444 (e. g. iFinke et al.l toid ): lOominguez et al.l J2OI1I ): lAleksic et aP feoilh l 



Notably, iGilmore et al.l ( 120121 ) shows that the intrinsic spectrum derived from deabsorption 
of 3C66A with the tentative redshift of 2; = 0.444 is the hardest of the deabsorbed VHE BL 
Lacertae objects. 

Figure 5 sho ws the VERITAS measured VHE spectrum of the blazar 3C 66A from 
Abdo et al.l ( 1201 ll ) (black solid line) when deabsorbed for the redshift upper and lower limits 
from this work. These deabsorbed spectra are calculated by multiplying the measured dif- 
ferential flux values by e^'-^'^-* for various EBL models. The resulting intrinsic flux estimates 
are then reflt with a differential power-law for the redshift lower limit (top) and 99% upper 
limit (bottom). The fltted intrinsic indices for both the lower and upper limits on redshift 
are su mmarized in Table 1. The hardest deabsorbed spectra result from the iFinke et al. 
( I2OIOI ) EBL model, but all fltted power-laws provide indices softer than the F = 1.5 limit 
(shown for reference in Figure 5 by the grey solid line at a comparable normalization to the 
deabsorbed spectra). The resulting indices are also below the Fermi-LAT measured index 
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ofr=i.8±o.i. 



180 4. Conclusion 

181 Observation of tlie z ~ Lya forest in tlie direction of the 3C66A with HST COS 

182 provides a direct lower and statistical upper redshift limit for the blazar. The detection of 

183 three clouds at z^bs ~ 0.3283, 0.3333, and 0.3347 provide the z = 0.3347 lower limit on 

184 the blazar redshift. Assuming that the incidence of Lyman absorption systems is a Poisson 

185 distribution in z, we can conclude that the blazar is likely to lie Zbiazar ^ 0.41 (99% confidence 

186 level) and exclude a z > 0.444 at 99.9%. 

187 Based on the assumption that the intrinsic index cannot be harder than F = 1.5, the 

188 redshift limits derived from the FUV observations do not place the blazar at a sufficient 

189 distance to utilize the observed VHE spectrum during an elevated state in October of 2009 

190 to constrain the EBL density. Moreover, the distance is not sufficient to extract an upper 

191 limit on the EBL density based on the similar assumption that the intrinsic VHE index is 

192 not harder than the Fermi observed index. 
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Fig. 1. — Keck/LRIS spectra (black) and error array (scaled by 50 x; grey) of the optical 
emission from 3C66A from September 2009 (top, relative high state) and October 2011 
(bottom, relative low state). The gaps in the spectra are due to the dichroic filter of the 
instrument. We have additionally cut the 2011 spectrum at 7000 A due to uncertainties 
introduced in calibration. All significant absorption features identified in the spectra are 
associated with the Earth or Milky Way. The details of thi s spec tral analysis for each of 
these observations are completed as described in lAbdo et al.l (1201 ll ). Even at this exquisite 
S/N (over 100 per pixel for both exposures) there are no features with which to place a 
constraint on the redshift of this blazar. 
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Fig. 2. — Detail of the COS spectrum of 3C66A in the regions where we identify Lya (top) 
and corresponding Ly/3 (bottom) absorption hnes for three gas clouds at ^abs ~ 0.3283, 
0.3333, and 0.3347. Absorption associated with Galactic Ni II is also labeled in the bottom 
panel. 

Table 1. Intrinsic indices (T) resulting from the deabsorption of the VERITAS observed 

spectrum reported in lAbdo et al.l (120111 ). Indices are calculated by taking the 
VERITAS-measured differential flux and flux errors and multiplying by e"^, where r is an 
energy and redshift dependent optical depth taken from the EBL models. The resulting 
flux in each bin is then fit with the differential power-law of the form dN/dE = {E/Eo)~^, 

where Eo is 250 GeV. 



EBL Deabsorbed Deabsorbed 

Model Index Index 

Used z = 0.3347 z = 0.41 



Gilmore et al. (2012'! 2.8±0.6 2.3±0.6 

Finke et al. (2010) 2.4±0.6 1.9±0.6 

Kneiske fc Dole (2010 2.8±0.6 2.4±0.6 

Dominguez et al. (2011) 2.6±0.6 2.1±0.6 
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Fig. 3. — Red portion of the G160M spectrum, redward of where we identify Lya hnes 
at Zabs ~ 0.33. All the labeled lines arise in the Milky Way. The lack of absorption of 
extragalactic origin places the redshift lower limit of 3C66A at -Zbiazar > 0.3347. 
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Fig. 4. — Top Distribution of the number of lines detected in 1000 mock spectra for 0.335 < 
z < 0.444. Bottom The probabihty to observe no Lya hues if 3C 66A hes beyond given 
the expected number of Lya hues in the redshift interval z\\ < z < z^^ derived from Monte 
Carlo simulations (dotted line). 
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Fig. 5. — Deabsorbed spectra for 3C66A for the 2:11 of 0.3347 (top panel) and 99% confidence 
level Ziii of 0.41 (bottorn panel), where the observed VHE spectrum (black solid line) is taken 
from I Abdo et al.l (120111 ). with an index of F = 4.1 ±0.6stat for the applied differential power- 
law of the form dN/dE = {E/Eo)~^. For reference, a spectrum with an index of F=1.5 is 
shown as the theoretical limit for an intrinsic index, as explained in the text. The resulting 
indices for each redshift and model are summarized in Table 1. 



